screen=c(0,1,2,3,4)
inc.std=c(0.52, 0.46, 0.50, 0.56, 0.54)
spend.std=c(-0.04, -0.17, 0.21, 0.53, 0.54)
inc.spend=c(-0.05, -0.19, 0.21, 0.49, 0.53)

x.at=c(0,1,2,3,4)
y.at=c(-0.2, 0, 0.2, 0.4, 0.6)

png("figure4.png", width=5,height=5,units="in",res=1200)

plot(screen, inc.std, axes=F, type="n", ylab="Correlation between the two variables", xlab="Number of Screeners Answered Correctly", ylim=c(-0.2, 0.6), xlim=c(0,4))
lines(lowess(screen, inc.std), lty=6)
lines(lowess(screen, spend.std), lty=2)
lines(lowess(screen, inc.spend), lty=1)
axis(1, x.at)
axis(2, y.at)
legend("bottomright", legend=c("Corr Income Ineq and Std of Living", "Corr Spending and Std of Living", "Corr Income Ineq and Spending"), lty=c(6,2,1), cex=.66)

dev.off()

